Mathematical Medicine and Biology: A Journal of the IMA

Anion exchanger 4 (Ae4), encoded by solute carrier family 4 member 9 (Slc4a9), is required for normal salivary secretion, but how it supports fluid secretion remains poorly understood. Nearly a decade ago, we developed a mathematical model that attributed reduced secretion following Ae4 knockout primarily to Ae4 mediated monovalent cation transport. Here, we reassess that interpretation by developing a dynamic model of salivary fluid secretion informed by subsequent experimental findings. The model incorporates {beta}-adrenergic regulation of Ae4 and an additional apical chloride conductance shared by the control and knockout. To isolate the contribution of Ae4, the NKCC1 and Ae2 cycle rates calculated from the evolving control state are also prescribed in the knockout. Under these conditions, the knockout model secretes 20% less fluid over 600 seconds despite identical apical chloride conductances. The control begins with a larger intracellular chloride store and continues to acquire chloride through Ae4 during stimulation. The intracellular chloride difference increases during stimulation, whilst the weaker apical electrochemical driving force in the knockout reduces chloride export and fluid secretion. The analysis also identifies the alkalinity requirement of sustained Ae4 transport and quantifies the additional NKCC1 and Ae2 transport needed to compensate for its loss. These results support intracellular chloride availability as a mechanism linking Ae4 loss to reduced salivary secretion.

Buruli ulcer (BU) is a neglected tropical skin disease caused by infection with Mycobacterium ulcerans. The mode of its transmission into the skin has not been fully clarified, but it is believed to involve mechanical transfer, perhaps by insect bites. M. ulcerans secretes a toxin called mycolactone which causes necrosis of the skin and formation of large skin lesions. Understanding the development of these skin lesions over time is crucial to understand the pathogenesis of BU. We formulate a mathematical model to describe the space-time dynamics of M. ulcerans bacteria, mycolactone and skin cell densities as well as fibrin deposition. We investigate the relative contributions to ulcer formation of direct cytotoxicity of mycolactone and ischaemia following the loss of coagulation control. Our results are consistent with existing biological evidence that higher mycolactone concentrations accelerate cell death. Furthermore, the model indicates that the rate of wound expansion is primarily controlled by the diffusion and net birth rate of bacteria, which sets the overall scaling for the speed of wound expansion, along with two dimensionless constants that define the shape of the trajectory of wound edge position against time, namely the ratio of the mycolactone diffusion coefficient to the bacterial diffusion coefficient, and the ratio of the mycolactone degradation rate to the bacterial net birth rate. When mycolactone diffuses significantly faster than bacteria spread or when mycolactone degrades more slowly than bacteria reproduce, the initial development of the wound is faster. In contrast, the long-term wound expansion speed appears to be relatively unchanged across a wide range of the dimensionless parameters.

Fractional differential equations are increasingly used to represent temporal memory in tumour immune and pharmacokinetic models, but changing the derivative order without a consistent time scale can make comparisons difficult to interpret. We study a dimensionless, commensurate Caputo system for tumour, immune, normal cell, and drug states. We establish positivity, boundedness, the tumour invasion threshold, and local stability conditions. We then use a fractional Adams Bashforth Moulton predictor corrector method with an explicit mesh convergence study. Two conclusions correct potentially misleading interpretations of the model. First, the equilibrium locations are independent of the fractional order: decreasing the order changes transient time scales but does not by itself alter the model's asymptotic states. Second, after a unit bolus the drug equation has the exact solution C({tau})=E(-K{tau}). Its maximum is therefore the initial concentration, not a delayed peak. For 0<<1, the algebraic tail makes the infinite-horizon area under the concentration curve diverge. In the illustrative coupled simulations, the apparent late tumour reduction at small is consequently driven by persistent modelled drug exposure and must not be read as evidence of superior therapy. All parameters are dimensionless and illustrative; no patient data are fitted. Source code and generated tables are included to make every numerical statement reproducible.

Oncogene amplification on circular extrachromosomal DNA (ecDNA) has been linked to poor prognosis and higher treatment resistance in multiple types of human cancer. ecDNA are mobile genetic elements lacking centromeres that are partitioned unevenly into daughter cells at mitosis. While random segregation of ecDNA contributes to gene-copy-number heterogeneity among tumour cells, how ecDNA contribute to phenotypic heterogeneity, via mRNA and proteins, and how this affects targeted therapy outcomes is still not understood. In fact, cancer cell populations with ecDNA show remarkable heterogeneity in mRNA concentrations, which is not explained by gene-copy-number heterogeneity alone. In this paper, we use a stochastic model of gene transcription to mathematically model mRNA production and decay in cancer cells with ecDNA, extending a method that has been used to study phenotypic heterogeneity in eurokaryotic and prokaryotic cells. Using our model, we explore the space of possible mRNA distributions from a given ecDNA distribution, under different model parameter regimes. Our model is able to produce realistic mRNA distributions and we use it to estimate transcriptional parameters from single-cell measurements of ecDNA and mRNA counts in cell-line data. We discover a remarkable similarity in estimated parameters across cell lines and amplification types, indicating similar processes underlying transcription in ecDNA and chromosomal amplifications. We also model genotoxic therapy and use our model to predict the mRNA composition of persisting cells, showing the usefulness of such a modelling framework to understanding the origins of treatment resistance in ecDNA-containing cancers.

Host immunity acquired through prior infections is a key factor shaping the spread of emerging variants. However, accounting for complete infection histories in mathematical models is challenging because the number of possible histories grows exponentially with the number of variants. Here, we present a flexible framework for modeling the co-circulation of multiple variants emerging at arbitrary times while incorporating cross-immunity and variant-specific characteristics. We introduce three tractable approximations: truncation, which retains only infection histories up to a given number of strains; one variant at a time, which is exact when variants circulation does not overlap; and two variants at a time, which is exact when at most two strains co-circulate. We systematically evaluate these approximations across a range of cross-immunity levels and inter-variant emergence delays. Truncation approximations perform best under short emergence delays and strong cross-immunity, whereas the other approaches outperform them at longer delays. When applied to SARS-CoV-2 variant dynamics, the second-order truncation and the two variants at a time approximation accurately capture incidence patterns, whereas the other approaches perform well only for specific subsets of variants. These results identify the parameter regimes in which simplifying assumptions remain valid, provide guidance for model selection, and highlight the key role of infection history in shaping epidemic trajectories.

Allergic conjunctivitis is an inflammation of the conjunctiva caused by allergen; it is common disorder affecting up to 40% of the population. In this work, we study seasonal allergic conjunctivitis (SAC), also called "hay fever eyes", which is caused by exposure to airborne pollens. We develop a mathematical model quantifying the ocular immune system response to the allergens. First, we present a simplified qualitative description of the immunopathogenesis of SAC. Then, we express each chosen immunopathological mechanism mathematically to construct a system of thirty-one ordinary differential equations. We compare summary statistics of the predicted observable immune signals to experimental measurements and find our model captures key qualitative features of SAC progression. We then compare our predicted time series of histamine concentration to symptom scores and find a strong correlation suggesting the model predicts relevant clinically trends. Next, we calibrate the model through multi-step process. We find the most influential parameters are the production and depletion rates of IL-4, and the production rates of IL-5 and IL-8. These cytokines are targeted in treatments for asthma, atopic dermatitis, and severe eosinophilic associated disorder and suggest potential therapeutic targets for SAC. Our calibrated model mimics most of the summary statistics of the experimentally observable immune signals with discrepancies for IL-5 and IL-13 indicating that additional immunopathological mechanisms could be important.

Background. Classical macroparasite models commonly represent worm burden heterogeneity through prescribed Poisson or negative-binomial (NB) distributions, while transmission dynamics are formulated in terms of aggregate quantities such as mean worm burden and prevalence. We revisit these assumptions using a stratified worm burden (SWB) framework that explicitly tracks the distribution of parasite burdens across heterogeneous host classes and allows demographic turnover, transmission, and interventions to be analyzed within a unified model. Methods. We derive equilibrium distributions of the SWB system under homogeneous and heterogeneous infection processes and compare them with the classical Poisson and NB distributions used in macroparasite epidemiology and control policymaking. The framework is then extended to include host demography, coupled human-snail transmission dynamics, and periodic mass drug administration (MDA). Corresponding SWB and NB-type MacDonald models are calibrated to comparable endemic prevalence and mean worm burden, allowing direct comparison of potential transmission breakpoints and intervention outcomes. Results. In the absence of demographic turnover, the SWB framework naturally generates Poisson and NB equilibrium distributions as limiting cases of homogeneous and heterogeneous infection processes. Incorporating demographic recruitment and mortality alters these benchmark distributions, producing endemic burden structures that differ systematically from their classical counterparts. Analysis of coupled human-snail systems reveals that transmission breakpoints and bistable equilibria occur only within a relatively restricted region of SWB parameter space. In contrast, comparable NB-MacDonald formulations exhibit breakpoint behavior over substantially broader parameter ranges with the net effect of overestimating the likelihood of interruption of transmission. Although calibrated models can reproduce near identical endemic prevalence and mean worm burden, they often generate markedly different responses to MDA, including differences in rebound dynamics, persistence, and elimination thresholds. Conclusions. The SWB framework links classical statistical descriptions of worm aggregation to mechanistic burden-resolved population dynamics, clarifying how endemic distributions emerge and how they are modified by demography and control interventions. Our results suggest that transmission breakpoints are not a universal feature of helminth systems but depend strongly on model structure, and that predictions of elimination and post-treatment recovery may be considerably more model-dependent than implied by prevalence-based calibration alone. The SWB approach therefore provides a flexible alternative framework for evaluating parasite transmission dynamics and predicting the outcomes of control strategies.

The aim of this study was to develop a framework for assessing the incremental value of decisions on drug amount, timing, input structure, cessation and restart. The LOGOS (Longitudinal Optimization of Global Outcomes and Strategies) approach combines matched nested regimen classes, mechanism-specific benchmarks and evidence requirements in a single longitudinal comparison. Thus, the contribution of additional regimen freedom is assessed separately from improvements due to different resources, constraints or evaluation horizons. Nested feasible classes have nondecreasing objective suprema, whereas a strict gain requires a feasible witness above a valid bound for the smaller class. The approach is illustrated through analytical examples of periodic smoothing, matched desensitization, relaxed exposure allocation, finite-course loading and safety after cessation. In a normalized finite Emax model, integrated benefit increased from 0.629606 for three equal boluses to 0.672945 after optimization of amounts, 0.677043 for an amount-and-time candidate and 0.679553 for a globally certified loading-maintenance-coast solution. The adjacent observed gains were 6.88%, 0.61% and 0.37%, respectively, and all schedules satisfied the full-tail injury bound. In addition, a periodic factorial example separated response curvature from a joint recovery and repair process, including a null case and results dependent on feasibility constraints. The LOGOS framework allows assessment of whether additional regimen flexibility improves benefit, provides no added value or restores feasibility. These dimensionless examples support conditional mathematical and model-based conclusions; clinical superiority and safe restart in patients have not been established.

We develop a general geometrical framework for growth processes in which the growth rate depends exponentially on the number of active growth-supporting components. The framework separates the fraction of active components from the size of the accessible component space, allowing qualitatively different growth laws to arise from the same underlying construction. Bacterial population growth provides an example of a fixed-size component space controlled primarily by activation and deactivation, whereas examples such as human population growth or speciation provide an expanding component space in which interactions between individuals promote further growth. This motivates a Replicative geometrical law, dX/dt= X exp({beta} X^{gamma}), whose exponent characterizes the scaling of the effective growth-supporting component space and is tested against historical human-population regimes and eukaryotic lineage-through-time data. For innovation indicators such as patents and arXiv manuscript submissions, where an existing item need not itself reproduce, the same construction motivates a Direct Geometrical production law without the prefactor X. The framework therefore provides a common description of fixed and expanding component spaces while distinguishing Replicative from Direct Geometrical growth.

Complex biological behaviours often emerge when a system is faced with mutually incompatible priorities. In such situations, the field of multi-objective optimisation can provide an informative and predictive theoretical foundation for biology. Here we adopt this paradigm in exploring the rich dynamic behaviour of mitochondria inside cells. Using plant cells as a model system, we use physical modelling to characterise the "morphospace" of possible mitochondrial behaviours, and single-cell microscopy with video analysis and network modelling to characterise collective mitochondrial dynamics. We show that wildtype Arabidopsis mitochondrial dynamics near-optimally resolve a tradeoff between maintaining physical spacing and supporting biomolecular exchange. With existing and new experimental data, we show that these dynamics adapt under mutational and chemical challenges to support a rebalanced, but still near-optimal, resolution to this tradeoff under different densities of the mitochondrial population. We also show how an assumption of multi-objective optimisation supports inference of biological mechanisms before any data are observed, and discuss the potential of this multi-objective optimisation paradigm to form a broader theoretical framework of spatial cell biology.

Devireddy et al. [1] previously showed that an enthalpy-based, macroscale model simulated tissue freezing histories are in close agreement with a coupled model incorporating cellular water transport and intracellular ice formation (IIF). That comparison, however, was performed for a limited set of tissue geometries, cooling conditions, and biophysical parameters. Here, we revisit that conclusion by reimplementing the published model from its governing equations and systematically expanding the parameter space to determine when microscale processes alter the predicted macroscale freezing response. We incorporate both surface catalyzed nucleation (SCN) and volume catalyzed nucleation (VCN) of intracellular ice with cellular osmotic dehydration and examine a broad range of tissue dimensions, cooling rates, convective boundary conditions, membrane permeabilities, cell sizes, activation energies, intracellular water fractions, and nucleation rates. Linearization of the cellular water transport equation about osmotic equilibrium yields an osmotic relaxation timethat can be compared with the local residence time of tissue undergoing phase change. Their ratio defines a dimensionless group, that organizes the transition between regimes in which the uncoupled enthalpy formulation and the coupled micro and macroscale formulation produce similar or substantially different results. Across the one-dimensional parameter space examined here, the two models remain within a prescribed thermal-error bound for the dimensionless group below approximately 0.1, whereas the discrepancy increases systematically above this threshold. Model divergence also depends on intracellular ice nucleation, i.e., when nucleation is sufficiently slow, cellular water transport can become rate limiting for latent heat release. Conversely, sufficiently rapid intracellular ice formation reduces the macroscopic discrepancy by providing an additional pathway for latent-heat release. Multidimensional simulations further show that geometry alone can move a nominally conventional freezing protocol into the coupled regime. In particular, multidirectional cooling can reduce the local residence time sufficiently to produce substantial differences in predicted frozen volume even when the corresponding thermal histories appear similar. Finally, we formulate an inverse approach for estimating IIF parameters from tissue scale thermal histories. Synthetic simulations indicate that parameter recovery is feasible within a finite range of cooling conditions, providing a testable strategy for future experimental measurements. Together, these results provide a mechanistic criterion for deciding when microscale biophysical processes must be retained in macroscale tissue-freezing models.

Cancer progression often requires the accumulation of multiple driver mutations, yet the same set of drivers may be acquired in different orders. How these alternative mutation-order pathways jointly shape tumor clonal structure remains unclear. We develop a multitype branching-process model in which malignant transformation requires two driver mutations and distinguish malignant cells both by their mutation order and by the independent transformation event that founded their clone. Under a successive exponential approximation, we establish point-process limits for the pathway-specific clone-size processes and derive a closed-form expression for the limiting expected Simpson's index of the combined malignant population. When the two mutation orders produce malignant cells with the same net growth rate, the limiting index separates into effective pathway weights, determined by mutation rates and birth-death dynamics at preceding stages, and within-pathway concentration terms, determined by intermediate-to-malignant growth-rate ratios. This decomposition shows that a driver's effect on heterogeneity depends critically on when it is acquired. A strong driver acquired early expands the intermediate lineage and increases the supply of independent malignant founders, whereas the same driver acquired at the final transition strengthens the growth and age advantage of early-founded malignant clones. Under additive fitness effects, these counteracting mechanisms can produce a non-monotone relationship between selective advantage and clonal concentration. We further show that threshold-like non-additive fitness effects can generate highly concentrated malignant populations, while order-dependent terminal fitness leads the faster-growing pathway to dominate asymptotically. Together, these results reveal how mutation order, mutational accessibility, selection, and epistasis jointly determine lineage-level intratumor heterogeneity.

Trigger waves are positive feedback-mediated traveling waves that propagate without slowing down or diminishing in amplitude, enabling biochemical coordination over distances where diffusion alone would be insufficient. Apoptosis can spread as trigger waves of caspase activity, as demonstrated in Xenopus eggs and egg extracts. We measured apoptotic wave speeds in extracts at two cell cycle states and found that waves propagated more slowly in M-phase arrested cytostatic factor (CSF) extracts than in interphase extracts (21 versus 27 m min-1). Luther's equation, v=2{surd}D{kappa}, predicts wave speed from the diffusion coefficient D and the positive-feedback doubling rate {kappa}. However, D and {kappa} were comparable between the two extracts and could not explain the speed difference. This discrepancy highlighted a key limitation of the classical equation: it does not account for an activation threshold in trigger wave initiation, and in fact the threshold required to trigger apoptosis was found to differ between cell cycle states. We therefore extended Luther's equation to include such a threshold, obtaining v{propto}{surd}(AD/Uc ), where A is the reaction rate and Uc is the activation threshold. Perturbing apoptotic sensitivity confirmed the predicted threshold-speed relationship. These findings establish the activation threshold as a key determinant of trigger wave speed and provide a broadly applicable scaling law for predicting the spread of switch-like biochemical signals.

Nitrogen dioxide (NO2) is a widespread air pollutant associated with the risk and prevalence of various respiratory infectious diseases. Epidemiological studies have reported associations between NO2 exposure and COVID-19 incidence, and biological evidence indicates that NO2 exposure increases vulnerability to infection. Motivated by these findings and by the scarcity of mechanistic epidemic models that incorporate air pollutant data, here we introduce an NO2- and mobility-informed metapopulation model and study the first COVID-19 wave in Oklahoma. County-level NO2 concentrations enter the infection rates through a first-order parameterization, alongside high-resolution mobility patterns capturing heterogeneous inter-county connectivity. Using numerical simulations, we demonstrate that higher environmental sensitivity produces earlier and more synchronized outbreaks across counties, while lower sensitivity yields delayed, asynchronous trajectories. When fitted to county-level COVID-19 incidence data, our model outperformed the classic SIR model (fitted independently in each county) in over 75% of counties in Oklahoma across all standard error metrics. The improvements were even more pronounced in forecasting tasks across epidemic stages, with county-level win rates against the SIR model of 59.74% (surge), 87.01% (peak), 80.52% (post-peak), and 88.31% (decline) by nRMSE. To isolate the contribution of NO2, we compared our model against a mobility-only model with NO2 removed and homogeneous infection rates across counties. The NO2 term improved county-level forecasts during the surge (94.81% of counties by nRMSE), post-peak (88.31%), and peak (62.34%) stages, though not during the early and decline stages, suggesting that NO2 information can generate better forecasts than mobility alone during periods of high incidence. Our results demonstrate that environmental exposure can enter the transmission structure of a mechanistic epidemic model directly and that doing so yields measurable forecasting gains.

Personalized medicine aims to tailor treatments to individual patients, but its implementation is limited by scarce clinical resources and heterogeneous treatment effects. Electronic health records now offer longitudinal data that can support data-driven therapy selection. We focus on SGLT2 inhibitors, medications with broad metabolic and cardiovascular benefits, whose clinical impact varies significantly based on individual patient profiles. We propose a counterfactual Structural Causal Model (SCM) to estimate individual treatment effects (ITE) of SGLT2 inhibitors from observational data. The model assumes counterfactual stability and uses a Gumbel-Max sampling mechanism to capture stochastic treatment assignment while enabling differentiable inference. Optimal treatments are defined by maximizing the counterfactual value of future diagnostic outcomes. We evaluate the Counterfactual Gumbel-Max SCM using a doubly robust value estimator and show that, across diagnostic categories, it outperforms real-world prescriptions in predicted treatment value, highlighting subpopulations where SGLT2 inhibitors yield superior outcomes. By integrating causal inference with electronic health data, this framework provides individualized counterfactual reasoning about SGLT2 treatment effects and can support clinicians in identifying optimal therapies for type 2 diabetes and comorbid conditions.

Non-exponential residence times give rise to history-dependent or memory effects in epidemic dynamics, yet their mathematical formulations remain fragmented. Here we develop a unified framework for a history-dependent SEIR model with general latent- and infectious-period distributions. We first establish the equivalence between the stage-age and renewal formulations and show that both admit a unified transition-flow representation. We then introduce the concept of a transition-flow multiplier, which provides a direct link between residence-time distributions and transition flows. When the multiplier is a weighted sum of shifted power functions, the model admits an exact representation in terms of Riemann-Liouville or tempered Riemann-Liouville fractional operators, providing a first-principles foundation for fractional epidemic models. Finally, using Euler-Lotka equations and stochastic-order arguments, we characterize how the latent- and infectious-period distributions shape epidemic growth and decay. A longer latent period slows both early growth and late decay, while, at fixed means, greater variability in the latent period accelerates early growth but slows late decay, whereas greater variability in the infectious period slows both.

Understanding the relationship between weather conditions and infectious disease transmission has important public health implications; nevertheless, the functional form of this relationship often remains poorly understood. Humidity has previously been suggested to have a U-shaped effect on the transmission of respiratory pathogens, with optimal transmission occurring at both low and high levels. Here, to assess the impact of this functional form and of its potential misspecification, we first simulate transmission rates based on the U-shaped relationship and then test whether an exponential function accurately estimates these rates across a range of climates. We further use a transmission model to simulate case count data for a range of hypothetical pathogens, comparing outbreak dynamics under the U-shaped forcing function to those under the best-fit exponential approximation. In our simulations, the model fit for the exponential function is climate-dependent, performing best in very dry and very humid environments. In locations where specific humidity values centre around the vertex of the U-shaped function, the exponential model underestimates the amplitude and frequency of incidence peaks. These findings demonstrate the potential consequences of assuming an incorrect functional form for weather-infectious disease relationships and highlight the need for more careful consideration of this functional form in modelling practice.

Fibrotic strictures are among the most disabling complications of Crohn's disease, permanently narrowing the bowel and impairing motility, yet no approved therapy reverses them. Chronic inflammation alters pacemaker-network coupling, smooth-muscle excitability, and calcium-dependent contractility, while fibrosis thickens the bowel wall, narrows the lumen, and changes tissue mechanics. The relative contributions of these coupled electrical, contractile, and structural alterations to motility loss remain unclear. To address this gap, we develop an integrated electromechanical finite-element framework for fibrostenosing Crohn's disease that couples a fibrosis-driven growth model with a FitzHugh-Nagumo electromechanical model. A full-factorial 25 design of experiments is used to quantify the relative effects of electrical diffusivity, excitation threshold, peak active stress, wall stiffness, and hypertrophic remodeling on cyclic lumen-volume deformation. Motility is quantified by the standard deviation of lumen volume over one contraction cycle. Within the parameter ranges examined, increased wall stiffness emerged as the dominant contributor to motility loss, followed by impaired smooth-muscle contractility. Changes in excitation threshold, hypertrophic remodeling, and electrical diffusivity produced substantially smaller effects. Pairwise interactions were small relative to the dominant main effects, indicating that the mechanisms contributed largely through their individual effects. Our findings suggest that limiting wall stiffening while preserving smooth-muscle contractile function may provide a therapeutic strategy for maintaining intestinal motility in fibrostenosing Crohn's disease.

Mechanistic models of infectious disease transmission provide a framework for understanding and predicting disease dynamics, but statistical inference is often computationally prohibitive. This challenge is especially pronounced for models with high-dimensional structures and many unknown parameters, where trajectory matching may require thousands of evaluations of an expensive likelihood function. We developed a surrogate-assisted optimization framework that uses Bayesian optimization to reduce the number of evaluations of the underlying mechanistic model required for parameter inference. We evaluated the approach using simulated fitting problems based on age-structured susceptible-exposed-infected-recovered models with 4, 8, and 16 age classes and a 210-compartment age-structured rotavirus transmission model. Applying surrogate-assisted optimization substantially reduced the number of likelihood evaluations required to identify good parameter estimates, with larger improvements as model dimensionality increased. In the rotavirus model, the surrogate Bayesian optimization approach reached the benchmark after approximately 0.7 hours, compared with approximately 9.4 hours using traditional Nelder-Mead optimization. These results demonstrate that surrogate-assisted optimization can substantially reduce the computational burden of trajectory matching for complex mechanistic transmission models, particularly for higher-dimensional inference problems.

We propose a stochastic epidemic model to describe the spread of avian influenza in a network of poultry farms. We consider three main paths of infection transmission between poultry farms: from nearby farms, mainly through airborne diffusion; from farms belonging to the same company, through shared veterinarians, forage providers and similar; from undetected small farms or wildlife. From the parameter estimation, based on a modified Expectation-Maximization algorithm, we infer that approximately 63% of the farms were infected from nearby ones, with an infection force declining with distance; of these, more than one third belonged to the same company of the estimated infector. About 20% were infected from premises belonging to the same company but farther away than the distance threshold of 2 km; the remaining ones from wildlife or unidentified sources. These estimates have been validated by a comparison with genetic data, available for a subset of the farms: the genetic distance between two farms identified, with high probability, as an infector-infectee pair is much lower than between random pairs. These results may help in the implementation of tailored prevention and control measures for future outbreaks.